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Abstract. We show that the so-called Gayley-Owocki (Doppler) heating is important for the temperature structure of the wind 
of main sequence stars cooler than the spectral type 06. The formula for Gayley-Owocki heating is derived directly from the 
Boltzmann equation as a direct consequence of the dependence of the driving force on the velocity gradient. Since Gayley- 
Owocki heating deposits heat directly to the absorbing ions, we also investigated the possibility that individual components 
of the radiatively driven stellar wind have different temperatures. This effect is negligible in the wind of O stars, whereas a 
significant temperature difference takes place in the winds of main sequence B stars for stars cooler than B2. Typical temperature 
difference between absorbing ions and other flow components for such stars is of the order 10'^ K. However, in the case when 
passive component falls back onto the star the absorbing component reaches temperatures of order 10® K, which allows for 
emission of X-rays. 

Moreover, we compare our computed terminal velocities with the observed ones. We found quite good agreement between 
predicted and observed terminal velocities. The systematic difference coming from the using of the so called "cooking formula" 
has been removed. 
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• ^ , 1. Introduction 

X 



Since the foundation of the theory of radiatively driven stel- 



lar wi nds by Lucy & Solomon (1970) and Castor, Abbott & 
Klein ( 1975 , hereafter CAK) many of the initial assumptions 
introduced by this authors were examined. To the most impor- 
tant on es bel ong the radial streaming approximation (Friend & 
Abbo tt [l98^ , Pauldrach e t al. |1986| ), the wind stability (Abbott 
1980|, Owocki & Rybicki|l984|), the limitations of the Sobolev 



approach (Poe et al. |l99Ci| , Owock i & Puis |1999| ), the thermal 
structure of the wind (Drew 1989) and many others. 

Another important assumption, studied already at the dawn 
of the radiatively driven stellar wind theory by Castor, Abbott 
& Klein (1976) is the condition of the one-component flow. 
They discussed encounters which transfer momentum received 
by absorbing ions (typically C, N, O, etc.) to passive, nonab- 
sorbing ions, mainly hydrogen and helium. They showed that 
for the high-density winds such encounters are not important 
for the overall dynamics of the wind and that high-density 
winds can be considered as one-component. However, for the 
low-density winds Springmann & Pauldrach (1992, hereafter 
SP) showed that momentum transfer between absorbing and 
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nonabsorbing plasma influences the wind thermal balance and 
even the wind dynamics. Thus, for the low-density winds the 
flow is essentially multicomponent. They proposed that the 
so-called "ion-runaway" may occur Based on the simplified 
theory of the multicomponent flow many of interesting results 
occurred. Porter & Drew (1995) re-examined model of wind- 



compressed disk in the presence of dynamical decoupling of 
absorbing ions and passive plasma. Porter & Skouza ( 1999| ) 
showed the possibility of formation of pulsating shells around 
stars with low-density radiatively driven wind, and Hunger & 
Groote (1999 1 explained the H/He abundance anomalies in Bp 
stars on the basis of helium decoupling. 

First detailed numerical models of multicomponent radia- 



tively driven stellar winds were presented by Babel (1995 
19961). However, Krticka & Kubat (|2000l hereafter KKO) 



showed that due to the functional dependence of the radiative 
force decoupling does not occur. Moreover, Krticka & Kubat 
(2001, hereafter KKI) using nonisothermal multicomponent 



models concluded that winds of B stars are frictionally heated 
in such amount that the possibility of decoupling of absorbing 
ions from the passive plasma is excluded. 

The solar wind is well-known to posses large tempera- 
ture differences between electrons and protons. Such differ- 
ences were obtained also by Biirgi (1992), who used the three- 
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component models of the solar wind. So the natural question 
arises, whether similar temperature differences exist in the ra- 
diatively driven stellar wind or, in other words, whether the 
assumption of equal temperatures of all wind components is 
acceptable. In this paper we intend to answer this question. 

Any effect which deposits heat separately to individual 
component of the flow may influence our results. Thus, we 
shall include the effect of Doppler heating, introduced in the 
stellar wind domain by Gay ley & Owocki (1994, hereafter 
GO). Because it arises from the dependence of the radiative 
force on the velocity via the Doppler effect it deposits heat di- 
rectly to the absorbing ion component and thus, it can trigger 
the temperature difference between absorbing and passive ions. 

Proper treatment of ionization balance may be important 
for correct description of decoupling of individual components 
of the flow. Thus, we decided to compute electrical charges of 
individual components using adequate ionization balance for- 
mulas. 

2. Boltzmann equation with veiocity-dependent 
force 

The procedure of the derivation of the hydrodynamic equations 
from the Boltzmann equation for particle distribution function 
Fg of the particle s is thoroughly described in a number of text- 
books. However, it is commonly assumed that the Boltzmann 
equation can be written in the form (we use the Einstein sum- 
mation law) 



ot oxh 



fsh dFs _ 
ms dish V 



(dFs 



(1) 



coll 



i.e. it is assumed that the force f ^ is independent of velocity, 
which is acceptable when we consider gravitational or electri- 
cal forces. However, the latter assumption is not valid for the 
case of a wind driven by radiative force coming from line ab- 
sorption, which is strongly dependent on a velocity gradient. 
Therefore, we rederive here hydrodynamic equations without 
an assumption of the force independent of velocity. In this case, 
the Boltzmann equation for the one-particle distribution func- 
tion Fs of particles of type s is written as 



dF\ dF^ 
at oxh 



d ( fsh 



dis 



dF, 
dt 



—Fs - ^ . (2) 



coll 



Here ^sh (^ = 1, 2, 3) are the velocity components of individ- 
ual particles of a type s with mass nis, and fsh are the compo- 
nents of an external force acting on them. Right-hand side term 
expresses the effect of collisions. By a definition, the integral 
of the distribution function over the velocity space is a number 
density of s-particles. 



d^Fs 



(3) 



Now, the usual way to obtain hydrodynamic equations is the 
following. One multiplies the Boltzmann equation (^ by mul- 
tipliers nis, TTisish and ms£,sh£.sk and integrates it over the ve- 
locity space. For the discussion of the role of velocity depen- 
dent force in the Boltzmann equation we confine us to the left- 
hand side of the Boltzmann equation. The right-hand side (i.e. 



the collisional term) remains unaffected by a presence of such 
forces, so we assume in this section that the gas is collisionless, 
i.e. that the right hand side of the Boltzmann equation is zero. 

2. 1. Continuity equation 

Multiplying the Boltzmann equation (^ by and integrating 
over the velocity space we obtain the continuity equation 



— (nsms) + (nsmsVsk) = 0, 



(4) 



where Vsk are the components of the mean velocity of particles 
s, and 



UsmsVsk =ms d^^iskF, 



(5) 



Here, the force term disappeared through integration by parts 
and the collisional term is also zero. Apparently, velocity- 
dependent external force does not change the continuity equa- 
tion. 

2.2. l\/lomentum equation 

Multiplying the Boltzmann equation (||) by rrisish and integrat- 
ing over the velocity space we obtain momentum equation 

9 , ^ d , , 

— (nsmsVsh) + t; — [nsirisVshVsk +Ps.hk) - 

ot OXk 



d^sfshFs = 0, (6) 



where Ps^hk are the components of the momentum transfer ten- 
sor, and 



UsmsVshVsk +Ps,hk = ms d^gishiskF, 



(7) 



Again, the velocity term was integrated by parts. Note, that due 
to the velocity dependence of the external force the term con- 
taining fsh cannot be moved in front of the integral. 

2.3. Energy equation 

Multiplying Boltzmann equation (^ by nisiskisk and integrat- 
ing over the velocity space we obtain energy equation 

a /I o 3 



dt\2 " 2^^ 

d fl n 3 1 

— ( -nsrusvsvs. + -PsVs. + VskPsM + ^Ps,kk, 

dUskfskFs - 0, (8) 



dx, \2 



where ps = i/'ips,kk is the scalar hydrostatic pressure and the 
last term was simplified using integration by parts. The equa- 
tion for the temperature can be derived by subtracting the mo- 
mentum equation (H) dot-multiplied by Vsh from the Eq. (M), 



d /3 



d 



dvsj 1 d 



+ Vsh f d^JshFs - [ d^sUfshFs - 0. (9) 
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Note that the last two terms cancel if the force does not depend 
on the velocity. For this "standard" case the external force does 
not give any contribution to net heating. 

2.3.1. Gayley-Owocki (Doppler) heating 

Let us explore now the last term on the left hand side of Eq. 
for the case of the force caused by absorption of radiation in 
spectral lines. This force depends on particle velocity through 
the velocity dependence of the line absorption coefficient ow- 



ing to the Doppler effect. Therefore, Gay ley & Owocki (1994 



hereafter GO) termed the heating effect by Doppler heating, 
but terming it Gayley-Owocki heating (or GO heating in the 
abbreviated form) might be more appropriate. Let us denote 
the heating term in the comoving fluid-frame as 



d$.iii,hfi,hF{. 



(10) 



Here the index i stands for "absorbing ions". Although this term 
is written in the comoving fluid-frame, the same expression 
holds in the non-relativistic case also in the observer frame. Let 
us assume complete redistribution and an angle independent 
opacity and emissivity in the atomic frame. However, the emis- 
sivity in the comoving fluid-frame is generally angle dependent 
due to the Doppler effect. Therefore, the radiative force acting 
on an atom in the comoving fluid-frame is 

h,h ^~ J dv <j> duj [K(n, v)I{n, v) - e(n, u)] nh, 

(11) 

where K{n, v) and e(n, v) are the absorption and emission co- 
efficients per unit mass, which can be expressed as 



(12a) 
(12b) 



K[n, v) — mp\v 4i • 

V c 



where (p{i') is the absorption (emission) profile in the atomic 
frame. After some rearrangement we obtain 



Qpo = ^ ^ diy (jiduj [Kl{n, v) - i 



1^0 



• (13) 



If for the calculation of the last integral we choose one of the 
velocity axes parallel to the direction of n, then the integrals 
over other two axes vanishes (we integrate odd function) and 
the Gayley-Owocki heating formula becomes 



^GO ^ ^ 
C 



J <j) dcj [Kl{n, ly) ~ e] X 

X I • (14) 



where is a velocity component in the direction of n. If we 
rewrite photon-line-of-sight velocity component as = wwth.i 



and use frequency displacement from the line center in Doppler 
units X — {ly ~ fo) /Aj^d, then the GO heating takes the form 

qGO _ / dx (p duj Kl(n, x) — Avjje x 

v'^^^-^dww Fi{vth.iw)'ip{w - x), (15) 



where we introduced the intensity /(n, x) as /(n, x) dx 
I{n^v) dv and the function 



is normalized according to Castor ( 1974 ) as 

'ip{x)dx = 1. 



(16) 



(17) 



Note that the thermal speed Vth.i is really ionic because it 
comes from the velocity distribution of absorbing ions (shall 
not be interchanged with vth, which comes from normalization 
of force multipliers). We neglect absorption in the resonance 
wings of the profile and approximated iIj{w ~ x) k S{w — x). 
Finally, we assume that the velocity distribution is given by the 
Maxwellian velocity distribution. 



Fi{vth,iw) 



Vth- 



The latter assumption was made purely due to simplicity rea- 
sons. Relaxing it could lead to interesting effects especially if 
the number of collisions i s not sufficient to ma intain an equi- 
librium (cf., e.g., Scudder 1994 , Cranmer 1998 ). However, we 
postpone the analysis of the non-Maxwellian effects to a future 
paper. 

Thus, the GO heating formula takes the form of 



-.GO _ TniUiVthA 



duj 



dxx(f>{x) X 
Kl{n, x) — Av^e 



where 



1 



(j){x) — —j=& 



(18) 



(19) 



Due to the symmetry of the absorption profile ( [l^ the product 
x(l){x) is an odd function and, thus, after integration over x van- 
ishes. This means that in the case of complete redistribution the 
process of emission gives no direct contribution to GO heating, 
i.e. 

Q^o = ^I^^I^^ ^ diu dxxcj){x)Ki{n,x). (20) 

In the static medium the product x(f){x)I{n, x) is an odd func- 
tion of X, thus, there is no GO heating effect in static stellar 
atmospheres. In the particular case of a spherically-symmetric 
stellar wind the expression for the Gayley-Owocki heating 
takes the form of 

GO ^ 27nw^ dxxq^{x)i{fi,x), (21) 

C J-l J-oo 

where fi — cos 6, which was actually used by GO. 
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2.3.2. Formula for Gayley-Owocki heating in tine 
stellar wind domain 

In the case of a two-level atom without continuum the solu- 
tion of the transfer equation in the Sobolev a pproxi mation is 
(Rybicki & Hummer |l978l Owocki & Rybicki |1985|, GO) 



(22) 



where Ic is the core intensity, D{fj.) is unity for > /i* and 

zero otherwise (/i* = (l — Rl/r^)^^^), core penetration and 
escape probabilities are given by 

1 



1 

2 



'-1 

respectively, and 



1 - e"^" 

T,, 



dx'(l){x'). 



(23) 
(24) 

(25) 



The Sobole v optic al depth is given by (Castor [1974] Rybicki 
& Hummer [l97^ 



where the variable <t was introduced by Castor (1974) 

dlnVri 



din? 



- 1. 



Inserting the solution of the transfer equation 
expression for the GO heating Eq.(El]) we obtain 



(26) 

(27) 
into the 



^GO _ 27rKpi?;tii,i/c 



dp 



D{p) 



/3 



X / da;x(/)(x)e-^^*(^). (28) 



The effect of line ensemble is usually described using the 
concept of a line-strength distribution function (CAK, Abbott 



1982, Puis etal. 20001 



dn 



dv 



10"cm-3 J V 
where normalization constant A^o is taken in the form of 

TVo = —{\-a)aa\-^. 



(29) 



(30) 



The GO heating formula for this line ensemble can be obtained 
by the integration of the heating term for one line (Eq, 28 ) over 
the CAK distribution function Eq.(^9[). In this case it takes the 
form of 



qGO _ Pi^^thWth.iiA^o ( Pel {Wmc) 



27rc22}ii?J V 10"cm-3 



dn K° 



dp 



D{p) 



dxa;(/)(x)e-^^*(^', (31) 



where Ic — iS.v-oL/ Att^R^ was used. Finally, applying substi- 
tution y — «;piVth7'/2)iWri preceding equation becomes 



-iGO _ PiVth,ikcFeL (1 - a) a / pc/ (Wuic) 



lOiicm-3 

2)iWri 



G{a,p,), (32) 



.CTePiVthr, 

where the function G{(j, p^) is given by the triple integration 



G{(J,p^) = / dyy" 
Jq 



dp 



dx x(f>{x) exp — 



1 + cr/i^ 



(33) 



and in the integrals for Pc and /3 Eqs. ( [23| , [23| ) the Sobolev depth 
shall be computed using 



(34) 



1 +(7^ 



instead of Eq. (pq). Contrary to the radiative force formula 
(Castor 1974) the GO heating formula depends on the absorp- 
tion profile. For the determination of the GO heating term we 
selected Gaussian profile (which comes from the Maxwellian 
velocity distribution). 



3. Model equations 

3. 1. Basic assumptions 

We assume that the stationary, spherically symmetric stellar 
wind consists of three components, namely absorbing ions, 
nonabsorbing hydrogen atoms and ions, and electrons, denoted 
by subscripts i, p, e, respectively. Each of them is described 
by a density pa, radial velocity Vra, temperature Tq, electri- 
cal charge qa — eza (where e is an elementary charge and 
Za denotes the ionization degree - may have a non-integer 
value), and particle mass nia- Subscript a stands for a = i, p, e. 
Contrary to our previous models (KKI), we allow for different 
temperature of each component and for radial changes of elec- 
trical charge. We assume that chemical composition is given by 
the factor z», which is a stellar metallicity relative to the solar 
value. 



3.2. Continuity equations 

In the case of a stationary spherically symmetric stellar wind 
each component is described by the continuity equation Eq.(0) 
in the form of 



J_d_ 

j,2 fj^j. 



(35a) 



where the term Sa accounts for radial change of mass-loss rate 
of individual components due to the ionization. Whereas for all 
types of ions the mass-loss rate is constant through the wind 
and thus number of these particles is conserved (Si — Sp = 0), 
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we account for the possibility of variation of electron num- 
ber via ionization and recombination. Because the total electric 
charge is conserved, 

^ dr \ ma J 

and continuity equation holds separately for aU components ex- 
cept electrons, we obtain electron continuity equation in the 
form of 



1 / 2 \ 



dZa 

^ — ' m„ ar 



(35b) 



Although inclusion of a term Sc into the electron continuity 
equation does not significantly alter the model, it is important 
to obtain well converged model. 

3.3. Momentum equations 

In the case of stationary spherically symmetric stellar wind the 
momentum equation Eq.(||) has the form of 



dVr 



rad 



dr 



9a 



— V KabG{Xab) 



Pa 



E 



qa_ 

ma 
Vrb 



E- 



\Vrb-Vra\'' 



(36) 



where square of isothermal sound speed is = kTa/ma, E 
is a charge separation electric field. Gravitational acceleration 
has the form g = GdJl/r^, where 9Jl is the stellar mass and G 
is the gravitational constant. The radiative acceleration acting 
on free electrons can be expressed as 



5c 



,rad 



(37) 



where the ratio of the radiative force caused by absorption of 
radiation by free electrons and gravitational force is 



A-KcGM 



L is the stellar luminosity and is the mass scattering coef- 
ficient of the free electrons (do not confuse it with Thomson 
scattering cross section). 

The radiative accelerati on acti ng on absorbing ions is taken 
in the form of Castor et al. ( |1975| ) 



rad 



fJeL ^ f Pel (Wnic 



2)i Airr^c 



2)i dVri 

aeVthPi dr 



(38) 



with force multipliers k, a, 6 after Abbott (1982). The finite 
disk correction factor (Friend & Abbott 1986, Pauldrach et al. 



1986) is 



/ 



(1 



, a+l 



(a + l)(l-Ai?)a(l + (7r 



(39) 



where W is a dilution factor. The thermal speed vth — 
-y/2fcTi/mp is, owing to the Abbott's normalization, hydrogen 



thermal speed. However, because vth physically describes ionic 
thermal speed, it depends on ionic temperature. Finally, 2)j is 
the photospheric ratio of the metallic ion density to the passive 
plasma density. We selected the same value as in KKI, namely 
2)j = 0.0127 (this value corresponds to the solar ratio of sum 
of densities of C, N, O, Fe to the density of bulk plasma). 

Constant of friction evaluated using Fokker-Planck approx- 
imation (cf. Burgers 1969) has the following form: 



Kab = nanb^^^ In A, 



(40) 



where and nb are number densities of individual compo- 
nents and mean temperature of both components 



Tab — 



TTlaTb + rUbTa 



TUa + mb 

The Coulomb logarithm is of the form 



In A = In 



ikT, 



kTr 



e \4TTne 



1/2- 



(41) 



(42) 



where n is the particle density (n = Up + ric + n-i). Finally, 
the Chandrasekhar function G{x), defined in terms of the error 
function erf (x) (Dreicer 1959[ ) is 



Gix) 



2x 



erf (x) — exp {—x^^^ 



The argument of the Chandrasekhar function is 

\Vrb - Vra\ 



•^ab 

where 



aab 



2 2k [niaTb + JUbTb) 

O^ab = ■ 

ma-mb 



3.4. Energy equation 



(43) 



(44) 



(45) 



Energy equation (^) in the case of a stationary, sphericall y sym- 
metric multicomponent flow has the form of (cf. Burgers 1969| ) 



-kVr 



' Tfla dr 



r dr 

1 ^ 2k {Tb - Ta) exp (-x^,^ 



-^Y^^ab 



b^a 



E 

b=^a 



nia + rub 

rrib 
rua + mb 



aab 



+ 



KabG{Xab)\Vrb-Vra\- (46) 



Two terms on the left-hand side stand for advection and adia- 
batic cooling, respectively. The right-hand side terms describe 
radiative heating/cooling, heat exchange by encounters of dif- 
ferent particles caused by unequal temperatures of the compo- 
nents, and frictional heating. 

There are two sources of radiative heating/cooling. First 
source are bound-free and free-free transitions and the second 
is Gayley-Owocki heating/cooling. 
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3.4.1 . Radiative energy term of bound-free and 
free-free transitions 

Bound-free and free-free transitions (which will be called 
"classical" radiative transitions) deposits energy directly to 
electrons. Therefore, this classical radiative energy term should 
be considered in electron energy equation. We decided to esti- 
mate the radiative heating/cooling term Q"^'^ using two mech- 
anisms only, hydrogen Lyman bound-free and free-free transi- 
tions. The detailed form of heating and cooling in the above 
mentioned transitions is nearly the same as in KKI and will 
not be repeated here (see also Kubat et al. 1999| ). The only dif- 
ference is that the temperature in these equations is now the 
electron temperature. at the base of the wind is taken as an 
emergent radiation from a spherically symmetric static hydro- 
gen m odel atmosphere for a corresponding stellar type (Kubat 
200l|). 



3.4.2. Gayley-Owocki lieating/cooling 

Contrary to bound-free and free-free transitions Gayley- 
Owocki heating/cooling deposits energy directly to absorbing 
ions. GO heating/cooling term has the following form (GO, 



-\rad 



GO 



GO 



PiVth.ikcreL (1 - a) a 



10"cm-3 



Gia,fi,), (47) 



where the function G{a-,fi^) is given by Eq.(|33[). Here ?;th,i 
is a thermal speed of driving ions, thus, vth.i = \/2kT]Jmi. 
Inspecting the GO heating term ( p7| ) we conclude that wind is 
heated via the Doppler effect by direct radiation whereas the 
Gayley-Owocki cooling is caused by the diffuse radiation. The 
sign of GO heating depends on the value of a variable a. In the 
case where cr = 0, when the expansion is locally isotropic, the 
GO heating and cooling vanishes. When the a is positive, the 
cooling term dominates and the wind is cooled by the GO heat- 
ing. In the opposite case, when the a is negative, the heating 
term dominates and thus, the wind is heated. 



3.5. Charge separation electric field 

The equation for charge separation electric field can be ob- 
tained directly from the third Maxwell equation, which in the 
case of spherical symmetry can be written as 



J_d_ 

J.2 fj^j. 



(48) 



Mihalas (1978, Eq. (5.46) therein), the ionization equilibrium 
in stellar winds can be approximated by 



1 ric Ua. 



iaj+1 2W Ua,j+i \2T:mkTj, 



3/2 



Tr 



1/2 



xe.p(|g), (49, 



where Xa.j is the ionization potential, Ua.j is the partition func- 
tion and Tr is the radiation temperature (we set Tr = |Teff). 
Partition function approximations are taken from Smith & 
Dworetsky (1988). Electrical charge of absorbing and passive 
ions are then derived using a formula 



Clearly, electron electrical charge is 



3.7. Critical points 



(50) 



1 everywhere. 



Critical points are points where derivatives of variables cannot 
be determined directly from differential equations. The deriva- 
tion of critical point conditions for our set of equations is sim- 
pler than that of KKI because we use different temperatures of 
each component and correct momentum equation for electrons. 
Due to these generalizations the set of critical point equations 
is not as coupled as it was in KKI and thus, the obtained critical 
point conditions are simpler. 

We write model equations in a simplified form, where we 
explicitly write only terms containing derivatives of individual 
variables and other terms are included into the terms Fi. Thus, 



the continuity equations (35a, 35b) are 



dr dr 



(51a) 



In the electron continuity equation (35b) we neglected the 
derivatives of ionic charge because their contribution to elec- 
tron continuity equation is only marginal. However, inclusion 
of such term influences critical point and regularity conditions 
for electrons only, which will not be used (see bellow). 

Similarly we can rewrite momentum equations (p6[). In the 
momentum equations of absorbing ions we shall linearize a 
term containing the velocity gradient. Note that because model 
equations are not quasi-linear (i.e. linear with respect to the 
derivatives of the independent variables), mathematically more 
correct method would employ some form of transformation to 
the quasi-linear form (cf. Courant & Hilbert 1962 ). However, 
because the results are essentially the same in this case, we 
present analysis of critical points in a simplified form. Thus, 
momentum equations are 



3.6. Ionic charge 

The ionization structure of stellar wind should be derive d usin g 
time consuming NLTE calculations (e.g. Pauldrach et al. 1994 ). 
Because we want to determine only a mean charge of selected 
elements, we use simpler approximate method. As described by 



dVr 



Vra- 



d9a 



ad 



dVr 



dr 



^ aidpa_ ^ o^dTa ^ 

d{dvra/dr) dr pa dr Ta dr 

= Fa,2{r,E,pb,Vrb,Tb,Zb). (51b) 



Similarly, due to the dependence of the Doppler term on the 
velocity gradient (in the Sobolev approximation) we shall write 
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energy equations in the form of 

3 al dTa 2 dvra 



Pa- 



dQa 



rad 



dVr 



d{dvra/dr) dr 

Fa,3{r, Pb,Vrb,Tb, Zb). (51c) 



The system of equations is closed by the equation for charge 
separation electric field, which has a simple form, 



dE 

— = Fi(r,pb,Zb). 



(51d) 



The system of equations ( |5T| ) can be simplified by insert- 
ing the derivatives of density from the Eq.(51a) and derivatives 
of temperature (51c) into the momentum Eq.(51b). We obtain 
modified linearized momentum equations 



However, as was discussed by KKI, this critical point condition 
is not met anywhere in the wind. Therefore, to obtain CAK type 
solution of the stellar wind we use another condition, which is 
connected to generalized Abbott waves (see KKI). Such condi- 
tion can be derived by multiplying the critical point conditions 
Eq.(|52[) by density of corresponding component and summing 
them 



EdVra 



dVr 



Pi 



' dr 



5ai_ 

3 Vra, 
Q - rad 



dgl 



,rad 



3uri d (dvr-Jdr) d {dvr-Jdr) 

= ^Fa{r,E,pb,Vrb,Tb,Zb). (56) 



rad 



5 a 



rad 



d{dVra/dr) iVra iv^ a d [dVr a / dr) 

dv, 



X 



dr 



Fair,E,pb,Vrb,Tb,Zb), (52) 



where Fa is a combination of Fa,i from Eqs. (51a), (51b) 
and (51c). Because the charge separation electric field equa- 



tion (51c) does not introduce any physically interesting critical 



point, the system of equations (^2|) consists of only three inde- 
pendent critical point equations. Each of them will be discussed 
separately in the following subsections. 

3.7.1 . Critical point of passive ions 

For the passive plasma the critical point condition ( ^2[ ) has a 
simple form 



2 ^2 



(53) 



In order to obtain continuous solution of model equations, at 
this point should the quantity Fa{r, E, pb,Vrb,Tb, Zb) vanish 
(regularity condition). This requirement is fulfilled if 



R-g+^E+-y^KpbG{x^b) 



3r 



2 1 1 



p6- 



\Vrb - Vrp\ 

2k in - Tp) {-4b) 



3 PpVrp ^ rup + mb apb 

- |— — E -T KpbG{Xpb)\Vrb - Vrp \ = 0. 

3 PpVrp ^ "^P + "^6 

(54) 

This equation is a generalization of well-known regularity con- 
dition for the coronal wind. 

3.7.2. Critical point of absorbing ions 

Critical point condition Eq.(|2|) for absorbing ions has form 



Vri 



dgl 



rad 



5 at 



dql 



rad 



d{dvri/dr) 3 Vri 3vri d (dvri/dr) 



= 0. (55) 



To assure that common equation of motion (obtained by sum- 
ming of individual momentum equations (^6|)) does not depend 
on the derivatives of any variables, a condition 



dVr 



dr 

dVri 

' dr 



5ai 

3 Vr, 



dql 



rad 



dgl 



rad 



3vri d (dvri/dr) d (dvri/dr) 



-0 (57) 



shall be fulfilled. This condition we use to fix the mass-loss 
rate. 

3.7.3. Critical point of electrons 

The last critical point condition for electrons has again very 
simple form 



(58) 



This critical point condition has similar form as the condition 
of nonabsorbing ions Eq.(^3|). Thus, regularity condition for 
electrons resembles the regularity condition for passive ions 
Eq.(i3), 



lOflp _j q, 
or mc p, 



E+ — Y,KcbG{x,b) 



Vrb ^re 



b^c 



\Vrb - Vre\ 



2 1 I 2fc(n-Te)exp 

o T= Z^^cb 

3 PcVrc ^ 



TOe + rUb Ueb 



I ^ — KcbG{Xcb)\Vrb - Vr 

3 PcVrc f-' rric + mb 

b^c 

2 1 



rad 



3 PcVr 



= 0. (59) 



However, numerical tests showed that electron regularity con- 
dition (|9|) is approximately fulfilled at the electron critical 
point if the zero current condition is used as a boundary con- 
dition. Thus, this condition should not be necessarily included 
into the set of model equations. 
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3.8. Boundary conditions 

3.8.1 . Boundary conditions for temperatures 

We assume that the flow at the inner boundary is in radiative 
equihbrium and that the boundary temperature of all compo- 
nents is the same, thus, we write boundary condition for tem- 
peratures in the form of 



Ta = Tc, a = i,p, 



(60a) 
(60b) 



Boundary values of ionic charges can be directly obtained from 
the condition of ionization equilibrium (49). 



Boundary value of ionic density is determined numerically 
to obtain CAK-type solution (see Sect. 3.9). Boundary electron 
density is calculated from the condition of quasi-neutrality 



(64) 



3.8.4. Boundary conditions for electric field 

Because we have not any critical point condition to determine 
the intensity of the electric field at the stellar surface, we used 
the condition of neutrality, which simply sets the gradient of 
the electric field at the stellar surface to zero (cf. Eq.Esl). 



3.8.2. Boundary condition for velocity 

Conditions ( jsj] , ^t] ) can be generally used to fix the boundary 
values of model quantities. However, inclusion of two inner 
conditions directly into model equations sometimes leads to 
numerical problems. Therefore, we use more secure method, 
which gives essentially the same results. 

We start to calculate our models at the passive plasma crit- 
ical point. Consequently, the boundary condition for the pas- 
sive plasma velocity is the critical point condition Eq.(53[). 
Boundary condition for the velocity of absorbing ions may be 
obtained from the passive plasma regularity condition Eq.(^4|). 
Because we suppose equal boundary temperatures of each 
component Eq.( |60a[ ), the regularity condition may be simpli- 
fied 



3.9. Numen'cal metliod 



lOai 
3r 



rr7__ n.. 



1 



2 1 ^ 

3 ppVrp 



nib 



pb) 



Vrb - Vrp 
\Vrb - Vrp\ 



6#P 



rrib 



-KpbG{Xpb)\Vrb ~Vrp\= 0. 



(61) 



The boundary value of electron velocity is chosen to fulfil 
the electron regularity condition Eq.(|9|) at the electron criti- 
cal point Eq. (p8|). As was already mentioned, this condition is 
approximately satisfied if the zero current condition 



'^7^ \ 



(62) 



is used as boundary condition for electron velocity. The latter 
condition (|2|) was applied in our models. 

3.8.3. Boundary conditions for density 

We write the boundary condition for the passive plasma density 
in the same form as in KKI, 



(63) 



Here we only newly introduced the relative abundance 
which accounts for different chemical composition in stellar 
atmospheres. 



We apply Henyey method (Henyey et al. 1964), which is a 
modification of the well-known Newton-Raphson method to 
solve equations described here together with the appropriate 
boundary conditions. We use essentially the same method as 
KKI, except that the vector of variables at each grid point d has 
the form of 



Ipd — {Pa,d, Vra,d, Ta^di Za,d, Ed, AVr^) j 



where the velocity difference 

AVr = Vr\ — Vr-r, 



e,i,p 



(65) 



(66) 



may be added to the set of variables to assure better conver- 
gence of the models. 

First of all we search for the boundary density po = pi ( i?* ) . 
We compute several wind models (each of them is a result of 
several Newton-Raphson iterative steps) for the region near the 
star for different values of po (for more details see KKO, KKI). 
We select such value of po which allows wind model to pass 
smoothly through the point defined by the Eq.(|57[) and to obtain 
CAK-type solution. After the appropriate value of po is chosen, 
we compute wind model downstream the point defined by the 
Eq.(p7|) again using several Newton-Raphson iterative steps. 

Detailed method of calculation of Gayley-Owocki heat- 
ing/cooling term is given in Appendix 

4. Results of calculations 



We computed several wind models for different stellar spectral 
types. Parameters of individual model stars are listed in Table |l]. 
Main sequence stellar parameters are taken from Harmanec 
(119881). For r Sco, e CMa, and /3 CMa we used the same 



parameters as SP and Cassinelli et al. (1995 1996), respec- 
tively. Note, that here we used slightly different parameters of 
T Sco than KKI. Force multipliers were adopted from Abbott 



(1982). The parameters of absorbing ions were selected in the 



following way. Because main sequence models were computed 
mainly for a demonstration of particular effects, it was suffi- 
cient to choose an ion which is simple enough and which de- 
scribes the basic line driving. So, we selected a carbon atom 
with rrii = 12TOp as a driving ion for them. On the other hand. 



J. Krticka, J. Kubat: Multicomponent radiatively driven stellar winds II. 



9 



3000 



2500 




100 



35000 
30000 
25000 
20000 
15000 
10000 





absorbing ions 




hydrogen 




electrons 




one-component 







10 

r/R. 



20 



50 



100 



Fig. 1. Upper panel: Comparison of the radial wind velocity 
of one-component (dotted line) and three-component radia- 
tively driven stellar wind models of an 06 star. Radiatively 
accelerated ions are denoted using dashed-dotted line, pas- 
sive plasma with full line and electrons with dashed line. 
Notice that all curves are very nearly the same. Lower panel: 
Comparison of temperature stratification of one-component 
and three-component models. Assignment of all curves is the 
same as for velocity. Notice that the curves describing temper- 
ature of individual components of a three-component model 
coincide, which is not true if we compare one and three com- 
ponent models. 



a driving ion of individual giant models was selected more 
carefully with respect to the stellar type. For e and (3 CMa we 
selected iron as a driving ion. The effective temperature of r 
Sco is higher, thus, we again selected carbon as a driving ion 
for this star. We stress that the selection of driving ions does 
not influence the amount of radiative force. However, it affects 
the thermal balance of the wind (via the frictional and Gayley- 
Owocki heating). For comparison purposes we also computed 
nonisothermal one-component models (see KKI) of these stars' 
winds with the same stellar and wind parameters (however, 
without GO heating). 



3000 
2500 
2000 
1500 
1000 
500 





absorbing ions 
hydrogen 
electrons 
one-component 




100 



Fig. 2. The same as Fig.|l] for an 08 star There is no signif- 
icant difference in the temperatures of each component. The 
wind temperature of the one-component and three-component 
models differ due to the Gayley-Owocki heating. 

4. 1. Winds with Gayley-Owocki lieating 

As was shown, e.g., by SP, high density winds are well cou- 
pled. For such winds the effect of frictional heating is negli- 
gible. Similarly, due to a high wind density the heat exchange 
between individual components is capable to maintain the same 
temperature for all components. However, for spectral types 
cooler than 06 a subtle effect of GO heating/cooling influences 
the temperature structure. This behaviour is displayed in Figsjl] 
and H for wind models of 06 and 08 stars. There is a large ve- 
locity gradient near the star, the variable a (Eq. |7|) is positive 
and thus, the wind is very slightly cooled by the GO cooling 
compared to the model without GO heating/cooling effects. On 
the other hand, in the outer parts of the wind the velocity gra- 
dient is lower, the variable a is negative and thus GO heating 
dominates (cf. Eqs. |3^, |3^. At the outermost parts of the wind 
the temperatures of models with and without GO heating are 
again nearly the same mainly due to the lowering of a stellar 
angular diameter 

As was discussed by GO, for stars with a lower density 
wind GO heating and GO cooling are much more evident. This 
can be seen in Figs.|], Q for the case of BO and B2 stars. As was 
shown by KKI, for such stars the frictional heating is negligible 
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Fig. 3. The same as Fig.|l| for a BO star. For this star Gayley- 
Owocki heating (in the outer parts of the wind) and cooling (in 
the inner parts of the wind) effects are important for temper- 
ature structure. Note that the frictional heating is negligible in 
this case and that the temperatures of particular components are 
nearly the same. 



and thus, changes in temperature stratification are caused only 
by GO heating/cooling. Due to the changed temperature strati- 
fication another effect becomes important. Because the radia- 
tive force in the CAK parametrization (Eq.^) depends explic- 
itly on the thermal velocity, higher wind temperature causes 
its lowering. Lower radiative force in the outer parts of the 
wind (above the critic al poin t) leads to lowering the outflow 
velocity (cf. Vink et al. 1999 ). However, the description of the 
dependence of the radiative force on the temperature via the 
Eq.(p^ is only approximate. Further calculations needed for 
better quantitative understanding the dependence of the radia- 
tive force on the temperature are currently under way and will 
be reported in future paper(s). 



4.2. Different temperatures of individual 
components 

For stars with lower wind density, individual wind components 
have different temperatures. This is shown in Figs.|] and ^ for 
B3 and B4 stars. 
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Fig. 4. The same as Fig|l| for a B2 star. Gayley-Owocki heating 
and cooling is greater than for a BO star - cf. Fig.|^. Frictional 
heating is negligible. The temperatures of particular compo- 
nents are nearly the same. 



Near the stellar surface the wind is relatively dense, the heat 
exchange between individual components is effective and tem- 
peratures of individual components are nearly equal. Compared 
to the model without GO cooling the wind temperature is 
slightly lower 

However, this is not the case in the outer parts of the wind. 
The wind is heated by frictional and GO heating there. Because 
the wind is more tenuous, the heat exchange between individ- 
ual components is not so effective as it is near the star and the 
temperatures of individual components differ. The temperature 
of absorbing ions is the highest and electron and passive com- 
ponent temperatures are nearly equal. 

There are three possible mechanisms which heat (or cool) 
individual components of the flow selectively and thus, which 
could make temperatures of particular components different. 
First, radiative heating/cooling (in our case made by bound-free 
and free-free transitions) deposits (or picks up) energy from the 
electron thermal pool. This causes electrons to incline to tem- 
peratures given by former models for whole wind - see Drew 
( 1989 ) for the one-component case and KKl when frictional 
heating is included. Second, the GO heating cools absorbing 
ions below the point where a = (see GO), and heats them 
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Fig. 5. The same as Fig.|l]for a B3 star. The wind is heated by 
both frictional and Gayley-Owocki heating in the outer parts of 
the wind. The temperatures of absorbing ions and electrons are 
nearly the same whereas the ionic temperature slightly differs 
mainly in the outer parts of the wind. 
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Fig. 6. The same as Fig.|l]for a B4 star The wind is heated by 
both frictional and Gayley-Owocki heating in the outer parts of 
the wind. Notice that the effect of heating is more pronounced 
than for the case of a B3 star (Fig. ||) and the temperature at its 
maximum is larger than for a cooler star The temperatures of 
absorbing ions and electrons are nearly the same. 



above this point. Finally, frictional heating itself deposits ther- 
mal energy unevenly. From the functional behaviour of the last 
right-hand side term of Eq.(^6|) we can infer that temperature 
increase is proportional to p^/ {rua + nii,). Thus, frictionally 
heated are mainly low-density components, i.e. electrons and 
absorbing ions, both via their collisions with passive ions. 

Another effect which influences the temperature balance 
is the heat exchange between components, described by the 
second right-hand side term of Eq.(p6|). Clearly, the heat ex- 
change depends mainly on the product of number densities of 
components. Thus, similarly to the differences in velocity, tem- 
perature differentiation takes place mainly for the low density 
wind. As discussed above, in such winds absorbing ions can 
be heated more than other two components. On the other hand, 
due to their large number densities, electrons and passive ions 
will share nearly the same temperature. All these effects influ- 
ence models in Fig.^ and|6| 

We determine the charge separation field directly from the 
Maxwell equation (psf). However, this has only marginal effect 
on the wind models, because our models tends to fulfil quasi- 



neutrality. 



The dependence of ionic charge for selected wind models 
is given in Fig.^ Clearly, because the ionization equilibrium 
depends mainly on radiative temperature, the ionic charge is 
nearly constant through the wind. 



Many B stars exh ibit U V-excess (eg. Cassinelli et al. |1995 



199q , Morales et al. |2001|). Note, that frictional heating and 



Gayley-Owocki heating of th e stel lar wind could be one of the 
possible explanations (Babel 1995). 



Frictional Ineating in one-component models: The fric- 
tional heating term in the energy equation of a one-component 
wind can be roughly estimated without using multicomponent 
models. In the ionic equation of motion ( |3^ all terms without 
radiative force and frictional term (corresponding to encoun- 
ters between passive component and ions) are small. We shall 
neglect them in this paragraph to derive an approximate expres- 
sion for frictional heating. The radiative acceleration term can 
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Fig. 7. The dependence of the ionic charge on radius for differ- 
ent wind models. Accelerated ions are denoted using dashed- 
dotted line. 

Table 1. Adopted parameters of model stars. 971 is the stellar 
mass in units of a solar mass, i?* is the stellar radius in units of 
solar radius. Toff is the star's effective temperature, k, a, and 5 
are radiative force multipliers, and is the metallicity. 



Stellar 


Stellar parameters 




Wind parameters 




type 


m 


R, 


Toff 


z. 


k 


a 


5 


mi 
mp 


(star) 




[Rq] 


[K] 












06 


31.65 


9.85 


41700 


1.0 


0.174 


0.606 


0.120 


12.0 


08 


21.66 


7.51 


35 600 


1.0 


0.166 


0.607 


0.120 


12.0 


BO 


14.57 


5.80 


29 900 


1.0 


0.156 


0.609 


0.120 


12.0 


B2 


8.62 


4.28 


23 100 


1.0 


0.377 


0.537 


0.091 


12.0 


B3 


6.07 


3.56 


19 100 


1.0 


0.477 


0.506 


0.089 


12.0 


B4 


5.12 


3.26 


17 200 


1.0 


0.365 


0.509 


0.105 


12.0 


B5 


4.36 


3.01 


15 500 


1.0 


0.235 


0.511 


0.12 


12.0 


r Sco 


19.60 


5.50 


33 000 


0.3 


0.113 


0.604 


0.095 


12.0 


eCMa 


15.2 


16.2 


21000 


0.18 


0.135 


0.561 


0.092 


55.8 


/3CMa 


15.5 


11.6 


23 250 


0.39 


0.125 


0.564 


0.099 


55.8 



be expresses as 
1 



9. 



rad 



(67) 



Here the subscript s denotes absorbing ions (note, that in these 
equations we allow for more than one type of absorbing ions). 
The frictional term in the energy equation of one-component 
wind can be approximated as 



absorbin, 
ions 

E 



Ps9s 



Vrp) , 



(68) 



absorbing 
ions 



where for the determination of the drift velocity from Eq.(|67|) 
an approximation of the Chandrasekhar function 



G(x) 



(69) 
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Fig. 8. The same as Fig.|l] for a B5 star The passive compo- 
nent decouples just above the stellar surface and subsequently 
falls back onto the star. Note that the temperature of the ionic 
component reaches the value of the order of 10® K. The tem- 
peratures of absorbing ions and electrons are nearly the same. 



can be used. We computed a one-component model for a B3 
star where the frictional term ( |68| ) was inserted into the en- 
ergy equation of one-component wind model (KKI, Eq.(40c)). 
Although the temperature was overestimated by 10% (com- 
pared to the correct three-component model) in the region 
where frictional heating is important, the difference of the ter- 
minal velocity was only marginal. Note, that the frictional heat- 
ing approximation Eq.(^8]) can be used even when there are 
more than one absorbing ion component. 

4.3. Backf ailing of hydrogen to the stellar surface 

In the case of the wind with the lowest densities the absorbing 
ions are not able to accelerate sufficiently the passive compo- 
nent of the wind. Thus, the passive component is not dragged 
out of the atmosphere and falls back onto the stellar surface (see 
Fig.^for a model of a B5 star wind). Such reaccretion should be 
studied using hydrodynamical calculations (Porter & Skouza 



1999). Moreover, Babel (1996 1 showed that the hydrostatic so- 



lution for passive plasma and the wind solution for absorb- 
ing ions exists. Probably, this type of solution is common for 
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Fig. 9. The same as Fig.|l]for a main sequence star t Sco. The 
wind is heated both by frictional and Gayley-Owocki heating. 



low-density winds, because Dworetsky & Budaj ( 200C| ) whilst 
studying Ne abundances in peculiar HgMn stars showed that 
in these stars the radiatively driven stellar wind with hydrogen 
mass loss rate larger than 10^^"* OJlo/yr is not present. 

Decoupling of velocities of absorbing and passive compo- 
nents is accompanied by decoupling of temperatures of these 
components (see Fig.||). This effect is caused by the depen- 
dence of the amount of heat transferred between individual 
components on the velocity difference (see Eq.(p^). Note, that 
the absorbing component attains temperatures sufficient to pro- 
duce X-rays. This effect can help to explain enhanced X-ray 
emission observed in mid- and late-B stars (Cohen et al. 1997[ ) 
which cannot be regarded as a consequence of standard radia- 
tion driven wind-shock mechanism. Another model for X-ray 
emission based on the shock decoupling was given by Porter & 
Drew (119951). 



4.4. Wind models of particular stars 

As was shown by KKl, the heating effect is pronounced in the 
wind of T Sco. Thus, we decided to recompute the wind model 
with the inclusion of GO heating. This model is shown in the 
Fig.^ Similar effects as in wind models of main sequence stars 
occur in wind models of giants. This can be seen in Figs. |lo| 
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Fig. 10. The same as Fig.[l| for a giant star e CMa. The wind 
is heated both by frictional and Gayley-Owocki heating. Note 
that the wind temperature at its maximum nearly reaches the 
effective temperature of the star. 

and |ll] for the wind models of (3 CMa and e CMa, respectively. 
For all these stars both frictional and GO heating are important 
for the temperature structure of the outer parts of the wind. 
For the star t Sco we used lower than observed value 



of metallicity. Contrary to Kilian (1994) who determined the 
value = 0.6 we reduced the metallicity to z^, ~ 0.3 to enable 
larger frictional heating. This change reflects mainly uncertain- 
ties of our model, because, e.g., our model with = 0.5 and 
metallic component described by iron ions instead carbon ions 
yields nearly the same velocity and temperature stratification. 

Unfortunately, existing measurements of the terminal ve- 
locity fo r this star do not allow to verify our models precisely. 
Abbott ( 1978 ) and Lamers & Rogerson ( 1978 ) determined 
Voo = 2000 kms^^, whereas Lamers et al. (1995) measured 
Voo — 1000 kms^^. However, all of them claim that their val- 
ues are uncertain. Larger values of Voc are supported also by a 
detailed UV-fit of Hamann (|1981|). 



For e CMa we used metallic ity = 0.18, a value esti- 
mated by Gies & Lambert (1992). Similarly to t Sco, available 
determinations of terminal velocity have lower quality. Abbott 
( |1978| ) determined Voo = 700 kms"^ and Lamers et al. ( |l995| ) 
measured = 600 km s^^. However, it is not clear whether 
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Fig. 11. The same as Fig.|l] for giant star (3 CMa. The wind is 
heated both by Motional and Gayley-Owocki heating. 



the apparent discrepancy of theoretical and observational ter- 
minal velocities is caused by the models or is due to the inac- 
curate measurements. 



According to Gies & Lambert (1992) we reduced the metal- 
licity of (3 CMa to the value — 0.39. To our knowledge, 
there is no measured terminal velocity for this star available 
in the literature. Note that for both e CMa and (3 CMa en- 
hanced wind temperature can help to explain observed UV- 
excess (Cassinelli et al. 



1995 19961. 



5. Comparison of terminal velocities 

In addition, we decided to compare our predicted terminal ve- 
locities with that measured by Lamers et al. (1995, hereafter 



LSL). They found discrepancy between theoretical values ob- 
tained from a "cooking formula" of Kudritzki et al. (1989 



hereafter KPPA) and their experimental values. We computed 
wind models of 06 - B5 stars for which LSL measured the 
terminal velocity. Parameters of each wind model are given in 
Tab. ^. Stellar parameters are taken from LSL, wind parame- 
ters are adopted from Abbott ( 1982 ). For many stars we found 
quite a good agreement between observed and predicted ter- 
minal velocities (see Fig.|l2| for comparison of predicted and 
observed values). Although some values of predicted terminal 



CO 




3000 



^oo (predicted) [km s" ] 



Fig. 12. Comparison of predicted and observational (taken 
from LSL) values of terminal velocities. Vertical lines denote 
uncertainty of observed values. Straight line is one-to-relation. 
For comparison, we plotted theoretical terminal velocities com- 
puted by LSL using "cooking formula" (crosses). 



velocities miss the measured value significantly (e.g. for the 
star HD 166596), it is evident that the overall agreement be- 
tween our predicted terminal velocities and the observed ones 
is much better than that of the "cooking formula" of KPPA and 
the systematic difference, which was previously attributed to 
an overestimation of a (by LSL), has been removed. 

However, there are still differences between observed and 
predicted values of Voo- There are three basic reasons for this 
discrepancy. First, rotation lowers the terminal velocity (cf. 
Friend & Abbott |1986|). However, Petrenz & Puis (|2000|) us- 



ing 2D models showed that the influence of the rotation on 
terminal velocity in many cases is only marginal. Second, our 
wind models (especially the radiative force) are constrained. 
Although we included physical processes that have not been 
included yet (frictional heating, Gayley-Owocki heating, mul- 
ticomponent nature of the wind), there are still limits. Our treat- 
ment of ionization is only approximate, the equilibrium is not 
determined consistently with radiation field. In addition, our 
models are not fully consistent with respect to the radiative 
force, a proper NLTE treatment of the radiative transfer prob- 
lem would be very useful. This two reasons causes that many 
of the terminal velocities are not within quoted uncertainties. 
However, we plan to improve our models in near future. 

Another source of differences may come from uncertainties 
of stellar parameters derived from observations. Note that, e.g., 
stars HD 106343, HD 148379, and HD 190603 have fairly the 
same parameters, however different observed terminal veloci- 
ties. 

The "cooking formula" of KPPA should be consistent with 
detailed calculations of Pauldrach et al. (1986) with an ac- 



curacy about 5%. However, our predicted termina l velo cities 
correspond to those computed by Pauldrach et al. ( 1986 ), too. 
Thus, there is not clear source of discrepancy between termi- 
nal velocities observed by LSL and predicted using formula of 
KPPA. We stress that the effect of frictional or Gayley-Owocki 
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Table 2. Stellar and wind parameters of 06 - B5 stars selected from LSL. 371 is computed mass-loss rate. Terminal velocities 
measured by LSL (column Voo (LSL)) are compared with theoretical values obtained by LSL using a "cooking formula" of 
KPPA (column Voo (KPPA)) and with predicted ones computed with an assumption of a nonisothermal wind model (colunm Voo 
(predicted)). 
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heating on the terminal velocity are negUgible for the models 
described in this section. 



6. Conclusions 

We computed non-isothermal three-component models of OB 
star winds with allowing for different temperatures of each 
component and with inclusion of the Gayley-Owocki (GO) 
heating/cooling. We showed that temperature differentiation 
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takes place in the winds of B stars starting from spectral type 
B3. The temperature of absorbing ions is of the order 10"^ K 
higher than temperature of other components whereas the tem- 
peratures of passive plasma and electrons is nearly equal. The 
main sources which trigger the temperature differentiation are 
GO, frictional, and radiative heating. 

Another important effect studied in this paper is the GO 
heating and cooling, which is important mainly for the low den- 
sity winds. We showed that this effect is a direct consequence of 
the dependence of the radiative force on the wind velocity. We 
derived the GO heating formula directly from the Boltzmann 
equation. More subtle GO cooling operates near the star at the 
wind base whereas the GO heating affects the flow mainly in 
outer parts of the wind. These effects become important start- 
ing from stellar type 06. Frictional and GO heating provides 
a possibility for an alternative explanation of UV-excess ob- 
served in some B stars. 

At the lowest densities either the passive component falls 
back onto the star or purely metallic wind exists. If the reaccre- 
tion takes place then ionic components is frictionally heated to 
the temperatures of orders millions K creating corona-like re- 
gion. This effect can explain enhanced X-ray activity in many 
of B stars. 

Finally, we compared our computed terminal velocities 
with that derived from observation. There is quite good agree- 
ment between them. The systematic difference between ob- 
served and predicted (by a "cooking formula" of KPPA) termi- 
nal velocities found by LSL was removed. However, we found 
no effects of frictional or GO heating in our sample. 
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Appendix A: Numerical caicuiation of G{a, /i*) 

The function G{a,ii^) is computed using numerical quadra- 
ture. Firstly, integral over x can be efficiently computed using 



a Hermite quadrature formula (cf. Ralston 1965). Quadrature 



weights and knots were computed using a subroutine IQPACK 
which is an implementation of method described by Kautsky 



& Elhay (1982). Satisfactory approximation can be obtained 



using 20 quadrature points. 

For an angle integration we used Legendre quadrature for- 
mula with 5 quadrature points. Again, quadrature weights and 
knots were computed using subroutine IQPACK (Kautsky & 
Elhay |1982[ ). 

Finally, the integration over y was performed using 
Simpson quadrature rule. The quadrature integral (0, oo) was 
approximated by (10"'^, 10^) and the Simpson integration was 
divided into subintervals of power 10, with 10 quadrature 
points in each of them. 



Numerical tests showed that temperature computed with 
described approximation of the GO heating has an error less 
than 1%. 
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